The magnetization process in the 2-dimensional Ji — J2 model 
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We study the a = J2/ Ji -dependence of the magnetization process in the Ji — J2 model on a 
square lattice with frustrating couplings J2 along the diagonals. Perturbation expansions around 
a = J2/J1 = and = yield an adequate description of the magnetization curve in the 
antiferromagnetic and coUinear antiferromagnetic phase, respectively. The transition from one phase 
to the other (0.5 < a < 0.7) leaves pronounced structures in the longitudinal and transverse structure 
factors at p = (tt, tt) and p = (0, n). 
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I. INTRODUCTION 

Magnetization processes in quasi one-dimensional 
quantum spin systems have been studied intensively dur- 
ing the last yearsL. The Oshikawa, Yamanaka and Affleck 
quantization ruleEl predicts possible plateaus in the mag- 
netization curve at certain rational values of the magneti- 
zation Af , which can be derived either from the geometry 
of the unit cell or a mapping onto a one-dimensional sys- 
tem with modulated couplingsB 

In the latter formulation, the plateaus appear at those 
values where the wave vector of the modulation coin- 
cides with the soft mode momenta of the unperturbed 
system. The existence of soft modes (zero energy excit^ 
tions) is guaranteed by the Lieb Schultz Mattis theoremcl 
for translation invariant, one-dimensional systems with 
finite-range couplings. A rigorous extension of this im- 
portant theorem to two and three-dimensional systems 
does not exist to date. For this reason, the Oshikawa, Ya- 
manaka and Affleck quantization rule cannot be applied 
in a straightforward way to predict possible plateaus in 
the magnetization curves of two- and three-dimensional 
systems. 

Experimental results, however, just concern com- 
pounds with a higher dimensional coupling structure as 
there are: 

(A) CsCuClj [Ref. ^ Here, a plateau has been found 
at M/Ms — 1/3 {Ms = 1/2 is the saturating magnetiza- 
tion). An explanation of this feature has been given in 
Ref. |. 

(B) NH4CUCI3 [Ref. 1^ Here, magnetization plateaus 
have been found for M /Ms = 1/4, 3/4. The compound is 
suggested to be built up as a two-dimensional structure 
of interacting two leg zig-zag ladders.!] 

(C) SrCu2{BO-i)2 [Ref. Magnetization plateaus 
have been observed for M/Ms = 1/3, 1/4, 1/8. The 
compound is modelled by the two-dimensional Shastry 
Sutherland model.EI M An explanation of the plateaus at 
M/Ms = 1/4, 1/3 has been discussed in Ref. ^ 

In this paper we are going to study the magnetization 
process in a two-dimensional spin-1/2 Heisenberg model 
with Hamiltonian 



with isotropic couplings 

i7,=^JKx,y)S(x).S(y), / = 1,2 



(1.2) 



for nearest and next-to-nearest neighbour sites x and y: 

■/l(x,y)=^y,x+i+'5y,x+2 (1-3) 
J2(X, y) = ^y,,+ i+2 + '5y,x+i-2 ; (1-4) 

1 and 2 denote lattice vectors in horizontal and vertical 
directions. For convenience we choose the nearest neigh- 
bour coupling to be one (Ji — 1, i.e. a = J2/J1 — J2) 
and use 



E 



e'P-5,(x), 



j = 1,2,3 



(1.5) 



for the Fourier transform of the spin operators on the 
square lattice. 

If we assu me periodic boundary conditions the Hamil- 
tonian (1.1) is translationally invariant. Moreover it 
commutes with the total spin operators S^(p = 0) and 
53(p = 0). We will classify the eigenstates |£',p, 3,33) 
by the eigenvalues: 

E oiHi+ aH2 , 

p = {pi,p2) of the momentum operator, 
s{s + 1) of the squared total spin S^(p = 0) , 
S3 of the 3-component ^3 (p = 0) . 

The magnetization curve M = M{B) is computed 
from the energy differences 

B{M = s/N, a) = E{s + 1, Ps+i, a) - E{s, p,, a) . 

(1.6) 

E{s,Ps,a) is the lowest eigenvalue of Hi + aH2 in the 
sector with total spin s. N is the total number of sites 
and Ps is the ground state momentum in this sector. 

Ps can be deduced from Marshall's sign rulellj in the 
limiting cases a = and a ^ 00. 

a) a = 



H = JiHi + J2H2 - BSsip = 0) 



(1.1) 



N 

Ps = (0, 0) if + * even 
N 

Ps = (tt, tt) if — + s is odd , 



(1.7) 
(1.8) 



1 



where the latter equations hold for clusters with size 
N = L X L, L even and periodic boundary conditions. 
Therefore, the transition between two subsequent ground 
s tate s as they enter on the left- and right-hand side of 
(1.6) is accompanied by a momentum transfer 



(7r,7r) 



6) a — > oo 

The Hamiltonian 



(1.9) 



(1.10) 



decays into two nearest-neighbour Hamiltonians on the 
even and odd sublattices, respectively. The nearest- 
neighbour couplings are defined here along the diagonals 
1-1-2, 1 — 2. Seen from the original lattice, the diago- 
nals are rotated by ±j. Moreover, the lattice constant 
increases by a factor \/2. 

We conclude from this that the momentum transfer 



between the two subsequent ground states in (1.6) is here 
q = Ps+i-p. = (vr,0), (0,^). (1.11) 



The change from (L9) to ( 1.11 ) signals a phase tran- 
sition fronipaptiferromagnetic to coUinear antiferromag- 
netic order .EJ 

Linear spin wave theory predicts that the regimes for 
antiferromagnetic and coUinear antiferromagnetic order 
are restricted to a < 0.4 and a > 0.55, respectively. 
In between., a phase with transverse disorder has been 
suggested.Ej 

We pursue the following strategy to study the impact 
of frustration on ground state energies and magnetization 
curves: We derive perturbation expansions in a and a~^, 
which are aimed to describe adequately the behaviour 
in the antiferromagnetic and coUinear antiferromagnetic 
phase, respectively. Comparison with computations on fi- 
nite clusters indicates that the perturbation results agree 
for a < 0.5 and a > 0.7, respectively. The regime inbe- 
tween which is not accessible by perturbation methods 
is of special interest. Here, we expecLpLhe emergence of 
plateaus in the magnetization curvcEjilj 

The outline of the paper is as follows: In Sec. II 
and Appendices ^ and ^ we present the results for the 
ground state energies obtained from perturbation expan- 
sions around a = and f3 — 1/a — 0. 

Free energies and magnetization curves are discussed 
in Sec. HI. Numerical results on the lowest frequency 
moments of the dynamical structure factor for momenta 
q = (tt, tt) and q = (tt, 0), (0, tt) are given in Sec. IV. 



II. FROM ANTIFERROMAGNETIC TO 
COLLINEAR ANTIFERROMAGNETIC ORDER 

In the antiferromagnetic phase a < ao(M) pertur- 
bation theory up to second order yields for the lowest 



eigenvalue of Hi + aH2 in the sector with magnetization 
M = s/N: 

E{s,ps,a) ^ N {ei{M) + ae2{M) + aH2{M)) ; 

a < ao{M) (2.1) 

where 



e,(Af) = l(0|i/,|0), J = 1,2 



(2.2) 



are the expectation values of Hi and H2 determined for 
the unfrustrated ground states |0) (i.e. H ~ Hi). 

The second order contribution 62{M) can be expressed 
in terms of the transition probabilities |(n|iJ2|0)p and 
energy differences i?„ — Eq between the ground state |0) 
and the excited states In): 



S2{M) 



-E 



En — En 



(2.3) 



The M-dependence of ei{M) and e2{M) is shown in 
Fig. |l](a) and |l|(b) respectively for system sizes = 4 x 4 
(with periodic boundary conditions) and A^ = 5x5=Fl = 
24, 26 (with helical boundary conditions). The data 
points for ei(M) nicely scale in M; deviations from scal- 
ing appear in £2{M) for smaller M values. 

Numerical results for the second order j-Gontribution 
S2{M) obtained with the recursion methodllZl are shown 
in Fig. ||. Scaling in M is realized for larger M-values 
(M > 0.3). For M = and M = 1/4, the TV = 4x4= 16 
data significantly deviate from the larger system results 
with TV = 5 X 5 =F 1. We suggest that the deviations from 
scaling in e2(M) [Fig. |l|(b)] and fe/M) [Fig. H arise from 
peculiarities of the 4x4 systemllj and that the thermo- 
dynamical limit is fairly well approximated by the larger 
system results. 
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FIG. 1. M-dependence of ei,2(M") for system sizes 
= 4 X 4 (square lattice) and N — ±1 (helical boundary 
conditions) 
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FIG. 2. M- and A'^-dependence of the second order contri- 
bution (52 (M) 

In the coUinear antiferromagnetic phase a > ao{M) 
perturbation theory in (3 — 1/a to second order yields 
for the lowest eigenvalue of Hi + aH2 in the sector with 
magnetization M = s/N: 



E{s,ps,a) = N {2M^ + aei{M) + a-^5-i{M)) , 



a > aa{M) . 



(2.4) 



The perturbative results to first orde r - i.e. the first 
two terms on the right-hand side of ( |2.4[ ) - are derived in 
App. A. Note that this result can be expressed in terms 
of the ground state energy per site ei(M) of the nearest- 
neighbour model H = Hi ! The second order contribution 
(5_i(M) is computed in App. B. 

We would like to stress at this point, that the pertur- 
bation expansions (2.1) and (2.4) hold in the thermody- 
namical limit. 

The determination of the coefficients ei(M), £2(^^)7 
(52 (M) and S^i{M) on finite clusters might be affected by 
finite-size effects; their magnitude is illustrated in Figs. 
|l|, I and |l^ and can be reduced by a computation on 
larger clusters and by a systematic finite-size analysis. 
The coefficients ei{M) and e2(M) are related to spin-spin 
correlators over nearest and next-nearest neighbour sites 
in the unfrustrated system (a = 0). Here, the finite-size 
effects are small and well under control. 

In Fig. ||(a) we show the a-dependence of ground state 
energies per site e{a,M = 1/ 4) as t hey follow from the 
perturbation expansions (2^) and (2^). The dashed 
lines represent the contributions to first order, which is 
fixed by the coefficients 



eiiM = 
62 



1/4) 

-5-1 



-0.34, e2(M = 1/4) = 0.46 



0. 



These values can be read off Figs. ||(a), (b). The sohd 
lines include second order contributions. Indeed, second 
order contributions seem to be relevant only in the an- 
tiferromagnetic phase (2.1). They are negligible in the 
collinear antiferromagnetic phase (2.4). The solid dots 
represent numerical results for a 4 x 4 system. Compar- 
ison with the perturbative results reveal small finite-size 



effects in the antiferromagnetic phase ( |2.lD , however large 
finite-size effects in the collinear antiferromagnetic phase 
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FIG. 3. Comparison of the a-dependence of the ground 
state energy per site e{a, M) - M = 1 /4 (a ), 5/12 (b) - and 
the perturbation theory results (^.l|), (2.4). The numerical 
data for M — 1/4 was taken from a 4 x 4 system and for 
M — 5/12 from a 6 X 6 system, both with periodic boundary 
conditions. 



The comparison of the perturbation expansions (2.1) 
and (2.4) with numerical results for M — 5/12 on a larger 
system with 6 x 6 = 36 sites is shown in Fig. ||(b). Note 
in particular that the deviations in the collinear antifer- 
romagnetic phase are smaller here, owing to the larger 
system size. In this case, the coefficients turn out to be 

ei{M = 5/12) = 0.17, e2{M = 5/12) 0.49 



The first order perturbation expansion ( ^.4| ) in the 
collinear antiferromagnetic phase a > ao{AI) predicts 
the vanishing of the a-dependence for AIq ~ 0.37, where 
ei(Mo~0.37) = [cf. Fig. |l](a)]. 



e(Mo, a) = 2Mq for a > aoiMo) . 



(2.5) 



Fig. |^(a) shows the Af -dependence of e(M, o;) on a 6 x 6 
system for various values of a. All curves with a > 0.6 
- where the system is in the collinear antiferromagnetic 
phase - coincide at Mq ~ 0.37. Data points for a = 
(open circles) - i.e. in the antiferromagnetic phase - do 
not share this property! 

The second derivative of e(M, a) with respect to M, 
computed from finite system results with = 36, AM — 
1/N: 



3 



c(2) 



9M2 



^(e(Af + AM,«) + 

+e(M - AM, a) 2e(M, a) ) (2.6) 



is shown in Fig. ^(b). 

The data points with M = 1/2 - = 36 (one 

spin flipped) have a dip at a = 1/2 indicating that the 
second derivative might vanish here in the thermodynam- 
ical limit N ^ oo. Indeed, we see from a Taylor expan- 
sion around M — Mq: 

e(M, a) = e(Afo, «) + e^^' (M - Mo) + ^e^^HM - M^f 



^i.e(3)(M-Mo)^ + 



(2.7) 
(2.8) 



M=Mo 



that the vanishing of the second derivative e^^^ — 
induces a square root singularity in the magnetization 
curve 



\M - M, 



01 



.(3) 



^\B-Bo\'/\ 



(2.9) 



provided that the third derivative e^^-* does not vanish. 
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FIG. 4. Ground state energy per site - e(M, a) vs. M - 
(a) and its second derivative - d^e/dM^ vs. a - (b). 

Under this assumption the numerical results for the 
magnetization curve at a = 1/2 and M ^ 1/2 on system 
sizes 5 X 5 ± 1 = 26, 24, 7 X 7 ± 1 = 50, 48 with hehcal 
boundary conditions could be described in Ref. |l^ by the 
square root (^). j— , 

However, it has been pointed out by HoneckerE3 in a 
recent publication, that numerical results on systems of 



size 6x6, 8x8 with periodic boundary conditions indicate 
that the magnetic fields near saturation 



B{M = n/N), a ^ 1/2, iV) = 4 

for n = 0,l,2,..L, N = 



(2.10) 



are independent of n. This means that the system jumps 
discontinuously from the magnetization M — 1/2 - (1 — 
L/N) into the fully magnetized state M = 1/2. The size 
of the jump AM = L/N vanishes in the thermodynami- 
cal limit, so one expects a smooth approa ch to saturation 
in this limit. Moreover the behaviour (2.1C) suggests, 
that all the derivatives (2_^) near saturation Ms — 1/2 



vanish, signalling an essential singularity. 

Going to a smaller value of M (M = 14/36 = 0.39) 
we find again a dip in the second derivative. However, it 
is very difficult to locate the dip position with our finite 
system results. 



III. FREE ENERGIES AND MAGNETIZATION 
CURVES 



According to eq. (1^) the magnetization curves follow 
from a minimum of the "generalized free energy" per site 



f{M, a) = e(M, a) ~ M ■ B 



(3.1) 



with r esp ect to M at fixed B. The perturbation expan- 
sions (2J) and (2^) therefore yield for the magnetization 
curves in the antifcrromagnetic and coUinear antifcrro- 
magnetic phase: 



BiM,a) = ha ha 

^ ' ' dM dM dM 



B{M, a) = AM 



dei 
'dM 



dM 



(3.2) 



(3.3) 



In Fig, [sl we present on the left-ha nd sid e th e free 



energies (3.1) - as they follow from (2.1), ( |2.4| ) - for 
a = 0.5, 0.6, 0.7. 

For a = 0.5 [Fig. |(a)] we are definitely in the anti- 
ferromagnetic phase, since the corresponding free energy 
here (solid line) is below the one in the coUinear antifcr- 
romagnetic phase (dashed line). In order to demonstrate 
the quality of the perturbation expansion (2.1), we also 
plotted the first order contribution alone (dotted line) 
and the numerical data (open circles) on a finite system 
with iV = 5x5 — 1 = 24 sites (with helical boundary 
conditions). Good agreement is found for M < 0.3. For 
larger M-values, the numerical data are below the per- 
turbative results. 

The corresponding magnetization curves at a = 0.5 
can be seen on the right-hand side. Again, the first per- 
turbation expansion ( |3.2| ) (solid line) follows the numer- 
ical results (step function), whereas the second pertur- 
bation expansion ( |3.3D (dashed curve) lies above the nu- 
merical results. 
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FIG. 5. Free energies and magnetization profiles for 
a = 0.5(a), 0.6(6) and 0.7(c). Solid and dashed curv es s liow 
tfie second order perturbation results [cf. (3.2) and (3.3)] in 
the antiferromagnetic and collinear antiferromagnetic regime, 
respectively. The circles represent exact diagonalization re- 
sults onaAi' = 5x5 — 1 = 24 lattice with helical boundary 
conditions. 

For a = 0.7 [Fig. |(c)] we are definitely in the collinear 
antiferromagnetic phase. The corresponding free energy 
(dashed line) is now below the one of the antiferromag- 
netic phase (solid line). The numerical data are close to 
the dashed line. In the magnetization curve (right-hand 
side) the perturbation expansions (^|^) , ( [3.3| ) and the nu- 
merical results (step function) coincide for M < 1/3. 

In between (0.5 < a < 0.7) we therefore expect the 
phase transition. In Fig. ^(b) we show the situation for 
a = 0.6. For M < 1/3 and M > 1/3, the free energies 
are minimal in the antiferromagnetic and collinear an- 
tiferromagnetic phase, respectively. The numerical data 
(open circles) on the finite system with N = 24 sites are 
lying below the perturbative results. The largest devia- 
tions are found in the collinear antiferromagnetic phase 
(A/ > 1/3) indicating that both perturbation expansions 
and finite system results become questionable. 
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FIG. 6. Behaviour of the magnetization curves of systems 
with periodic (A'^ = 5x5) and helical (A'^ = 5x5 — 1) boundary 
conditions around magnetization Af = 1/4,1/3 and frustra- 
tion a — 0.6 

For small values of M [M < 1/3) the numerical data 
for the magnetization curve (right-hand s ide) are well re- 
produced by the perturbation expansion ( |3.2| ) for the an- 
tiferromagnetic phase. At Af 1/3 a plateau-like struc- 
ture seems to evolve in the in the numerical data for 
N = 5x5— 1 = 24 with helical boundary conditions. 
Unfortunately, we cannot check whether this structure 
will survive on larger systems with helical boundary con- 
ditions like A^ = 7x7— 1 = 48. However, the structure 
seems to appear as well (at M = 0.34) on a square lattice 
iV = 5 X 5 = 25 with periodic boundary conditions as is 
demonstrated in Fig. pi 



IV. DYNAMICAL STRUCTURE FACTORS AND 
FREQUENCY MOMENT SUM RULES 

In this section wc will study the impact of the tran- 
sition from antiferromagnetic to collinear antiferromag- 
netic order on the dynamical structure factors: 

Sjj (q, u, M, a) = ^Yl -^(^^ -Eo-L^)\ {n\S, (q 



(4.1) 



and their frequency moments: 
Kl^^ (q. A//, a)= dco uj^'S,, (q, c^, Af, a) . (4.2) 



In particular, we expect signatures in the a- 
dependence of the structure factors at momenta q = 
(tt, tt) and q = (tt, 0). 

In Figs. 0and^ we show the a-dependence of the trans- 
verse and longitudinal static structure factors on a square 
lattice with = 6 x 6 = 36 sites: 

k'^^^ (q, M, a,N = 36) for M = — ■ ( 13, 14, 15, 16 

36 V 

(4.3) 
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FIG. 7. Zeroeth transverse frequency moment (transverse 
static structure factor) Kg''' {q, M,a) = 5'ii(q, M, a) for sys- 
tem size A'^ = 6 X 6 for different magnetizations and momenta 
q = (7r,7i')(a) and q= (0, tt), (tt, 0)(b). 

The transverse structure factors (j — 1) for the 
momenta q = (tTjTt) [Fig. 0(a)] and q = (tt, 0), 
(0, tt) [Fig. 0(b)] reveal a striking similarity in their 
a-dependence: All curves for q — (tt, tt) and q — 
(tt, 0), (0, tt) show a sensitive change in their M- and 
a-dependence when entering the regimes a > 0.6 and 
a < 0.5, respectively. 

In particular we see in our data, that 

(q = (tt, tt) , Af , a) ~ M/2 for a > 0.6 . (4.4) 

This behavi our is a consequence of the ground state 
wave function (A2). 

The longitudinal structure factor [(O) for j — 3] at 
the momenta q = (tt, tt) [Fig. ||(a)] and q — (tt, 0), (0, tt) 
[Fig. 11(b)] look similar in the small-a regime {a < 0.5): 
One observes a constant behaviour with a and a mono- 
tonic decrease with M. For large values a, the longi- 
tudinal structure factors behave quite differently: For 
q — (tt, tt) the curves with M = s/N, s even and s odd 
approach different limiting values, respectively. This fea- 
ture is not visible for q = (7r,0), (0,7r). 

Numerical data for the static structure factor on a 6 x 
6 system with, periodic boundary conditions have been 
presented inJlj At fixed magnetization M = 1/4 - where 
the authors see a plateau in the magnetization curve - 
the data points for the longitudinal structure factor at 
q — (tt, tt) and q ^ (tt, 0) show a broad maximum in 
the regime 0.55 < a < 0.67, which is also visible in our 
results for larger Af values [cf. Figs. ||(a), (b) in the 
interval 0.5 < a < 0.6]. 



FIG. 8. Zeroeth longitudinal frequency moment (longitudi- 
nal static structure factor) Kq^\ci, M,a) = Ssslq, M, a) for 
system size = 6 x 6 for different magnetizations and mo- 
menta q = (tt, 7r)(a) and q = (0, tt), (tt, 0)(b). 

Let us next turn to the first frequency roament. As was 
pointed out by Hohenberg and BrinkmanEil and MiillercHl 
it can be expressed in terms of ground state expectation 
values for two-fold commutators: 

Xp)(q, A/, «) = ^ E -^0- B{S3n - ^30)) 

71 

x\{n\SMm\' 
= l(0|5+(q)i/5,(q) - i/5+(q)5,(q)|0) 



N 
1 

2N 



(0|[[5+(q),i/],5,(q: 



(4.5) 



5*30 and S^n are the 3-component of the total spin in 
the ground state |0 > and the excited states \n >, re- 
spectively. The calculation of the commutators yields: 



-(1 - cosqi)Q(i) - (1 - cosg2)Q(2) 
~a (1 - cos((7i + (72)) Cj(i + 2) 
-a(l-cos(qi -q2))Q(i-2) 



MjB{M, a) 



(4.6) 



where 



X 

-(0|5,(x)5,(x + z>)|0)) 
for j = 1, 2, 3, i> = i, 2, i + 2, i - 2 and 
Af, = Af for j = 1,2, 



(4.7) 



Af, 



for j = 3 . 



6 



It should be noted that the dependence on the mo- 
mentum transfer q = (gi, 92) is complet ely g iven by the 
Fourier factors on the right-hand side of (4.6). The phys- 
ical meaning of the first frequ ency moment can be seen 
from the definitions (4.1) and (4.2). 




FIG. 9. Frequency expectation values {u))-i{q, M, a) for a 
= 6 X 6 spin system for different magnetizations and 
q = (tt, 7r)(a) and q — (0, tt), (tt, 0)(b). 



The ratio 



< LUj > (q, M, a) 



J = 1,3 (4.8) 



yields the average excitation energy in the spectrum 
of states which can be reached from the ground state 
|0) — \ps, Ms) by means of the operator Sj{q) j = 1,2,3. 
This ratio is shown for the longitudinal case j = 3 and 
magnetizations M ~ 1/36 • (13, 14, 15, 16) on a 6 x 6 sys- 
tem in Fig. ^. 

At momentum q = (tt, tt) [Fig. ^(a)], we observe a 
sharp drop in the width of the excitation spectrum at 
a « 0.5, i.e. close to the phase boundary. 

At q = (tTjO), (0,7r) [Fig. ^(b)] the longitudinal struc- 
ture factors for all M-values follow the same straight line, 
which ends at a minimum at a — 0.5. 



V. CONCLUSIONS AND DISCUSSION 

The magnetization process in the two-dimensional Ji — 
J2 model is well described in the antiferromagnetic phase 
{a < 0.5) by a perturbation expansion around a = 
and in the coUinear antiferromagnetic phase {a > 0.7) 
by a perturbation expansion around = 0. In these 
regimes, the magnetization curves are smooth, monoton- 
ically increasing and convex. The change in the magnetic 



order can be seen in the a-dependence of the transverse 
and longitudinal structure factors at momenta q — (tt, tt) 
and q — (tt, 0), (0, tt), respectively. 

Using frequency moment sum rules we have also stud- 
ied the variation of the average excitation energy (lo) 
(4.8) with a. Approaching the transition region from an- 
tiferromagnetic to coUinear antiferromagnetic order, the 
spectral weight is more and more concentrated around 
low frequencies. 

Between the antiferromagnetic and coUinear antiferro- 
magnetic phase - i.e. for 0.5 < a < Jd,7 - a phase of 
transverse disorder has been suggestedJla 

This regime is not accessible with the perturbative 
methods developed in this paper; nonperturbative effects 
- like plateaus - are expected here. 

To study these, we have to rely on numerical results on 
finite clusters. The results available so far reveal a strong 
dependence on the size of the clusters and the boundary 
conditions. For the moment we meet the following situ- 
ation: 

(A) At a = 1/2 the approach to saturation {Ms = 1/2) 
is no longer linear but appears to develop a singularity. 
The type of the singularity is not yet clear: 

On one hand, numerical data on square lattices 6x6 = 
36, 8 X 8 = 64 with periodic boundary conditions reveal 
the peculiar property (|2.10|) for the magnetic fields near 
saturation. If we ass um e in the thermodynamic limit a 
series expansion [cf. (2/7)] of the energies per site e(M, a) 
at the saturation point Ms = 1/2, infinitely many deriva- 
tives [cf. (2.^)] would vanish, indicating an essential sin- 
gularity. 

On the other hand, numerical results [Ref. |l^ on finite 
lattices 5 x 5 ± 1 = 26, 24, 7 x 7 ± 1 = 50, 48 with hcUcal 



boundary conditions do not share the property {2.10) for 
the magnetic fields near saturatiion M 1/2, a — 1/2. 
It was suggested in Ref. |l^ that the magnetization curve 
develops a square root singularity, wh ich follows if the 
second derivative e'^^^(Af, a = 1/2) in (2.7) vanishes for 
M 1/2, but not the higher ones e'^M,a = 1/2) 
fc = 3,4, .... 

(B) At a = 0.6 strong evidence for a plateau at 
M — 1/4 has been reported by Honecker et al.ll3G3 who 
performed a finite-size analysis of exact diagonalization 
results on 3 clusters 4 x 4, 4 x 6 and 6x6 with periodic 
boundary conditions. They do not find any signature for 
a plateau at M = 1/3. In our results obtained from a 
A^ = 5x5— 1 = 24 system with helical boundary condi- 
tions we do not find a plateau at M = 1/4; but there are 
indications for a plateau-like structure at M = 1/3. We 
are aware of the fact, that this discrepancy might come 
from special finite-size effects due to the helical boundary 
conditions, which only disappear for much larger systems 
(e.g. A^ = 7x7— 1= 48, ..). However, these are not yet 
accessible. 

In order to test the dependence on the boundary con- 
ditions we have compared the magnetization curves for 
A^ = 5x5 — 1 = 24 with helical boundary conditions and 
for = 5 X 5 = 25 with periodic boundary conditions. 
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The results (Fig. look similar. This would mean that 
the plateau structures in the magnetization curve has a 
sensitive dependence on the system size N ^ L x L; in 
particular results for L even and L odd might be quite 
different. Therefore, the precise determination of the 
positions and widths of the plateaus in the frustrated 
Heisenberg model is indeed a very delicate problem. 

A better theoretical understanding of the magnetic or- 
der and of the mechanisms which create the plateaus is 
really needed. 
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APPENDIX A: THE STRONG FRUSTRATION 
LIMIT 



As was pointed out in Sec. I [cf. ( |l.lO| )], the diagonal 
couplings build up two independent nearest-neighbour 
Hamiltonians on the even and odd sublattice, respec- 
tively. 

We start with the eigenvalue equations for these sub- 
lattices 

52 U(±)(a,p,)= aia + 1) U(±)(a,p,). (Al) 

53 } <^ } 

In particular we are interested in the ground state 
with definite (squared) total spin and total magneti- 
zation M = a/{N/2). In the thermodynamical limit 
N/2 — > 00 the|-ttround state momenta Po- follow Mar- 
shall's sign rulell^. If N/i -I- cr is even, the ground state 
momentum turns out to be Po- = (0,0); ii N/A + a is odd, 
one obtains Po- = (7r,7r). The lowest eigenstates of the 

Hamiltonian H2 = -ff j;"'''' -|- h[ with even total spin and 
zero total momentum p = (0, 0) can be constructed as 
a product of states of the lowest eigenstates (Al) of the 
subsystems: 

^{s,ps = (0,0), TV) = *(+)(s/2,p,/2,7V/2) 

xvI/(-)(s/2,p,/2,A^/2). (A2) 

The two subsystems defined on the even and odd sites 
can be transformed into each other via a translation of 
one lattice spacing. We can therefore impose the follow- 
ing identification on the wave functions with identical 
quantum numbers in ( [A2| ) 

*(-)(^,P^,^^/2) =Ti*W(a,p,,iV/2), a = 1/2. (A3) 



r,^'(s, p, = (0, 0),N) = *(s, p, = (0, 0),N) , (A4) 

i.e. they are momentum zero eigenstates. For the proof 
of this statement one has to realize that the eigenstates 
on the subsystems behave as follows under translations: 

Tir2^'(=^)(cr,p,A^/2) =e'P' • *(±'(CT,p,Ar/2) 
TiT2-i*(±)(ct,P,A^/2) = e'P' ■ ¥^^<7,p,N/2) (A5) 
T^¥^Ha,p,N/2) = e'(Pi+P2) . *(±)(a, p, iV/2) 



It should be noted that the product states (A2) where 
the total spin s is distributed to equal parts on the two 
sublattices indeed yields the ground state with energy 

Eis, p, = (0, 0),N) = 2E(s/2, Ps/2,N/2) . (A6) 

Any other distribution, e.g. with a total spin s/2 + As 
on the even and s/2 — As on the odd lattice will lead 
to an eigenstate with higher energy. This is a conse- 
quence of the fact that the lowest energy is monotonically 
increasing and a convex function of the magnetization 
M = s/N. 

The construction of eigenstates of i/^^' + h[ ^ with 
odd total spin s and total momentum p^ = (tt, 0), (0, tt) 
is more involved. In this case an equal distribution of the 
total spin on the two sublattices is impossible since the 
sublattice spin a has to be integer. The state with lowest 
energy is found with the ansatz: 

*(s,p, = (^,0),Ar) = 
^ (vE'W(a,p.,Ar/2)vI/(-)(a- l,p,_i,iV/2) 



*(+)(a-l,p._i,7V/2)«'(-)(a,p,,7V/2)) (A7) 



where 



(A8) 



Imposing the identification ( [A^ ) on the ground states 
with equal quantum numbers, one finds the following 
transformation properties of the states (A7): 



ri«'(s,p, = (^,0),7V) = -^'(s,p, - (7r,0),iV) 
r2*(s,p, = (^,0),Ar) = *(s,p, = (7r,0),iV) 



(A9) 



(AS) guarantees that the product states (A2) are in- 
variant under both [v = 1,2) translation operators: 



which means that the ground state momentum is Ps = 
(tt, 0). The degenerate ground state with momentum 
Ps) = (0,7r) is obtained if we substitute the translation 
operator Ti by T2 in (^). 

In the strong coupling limit a ^ 00 {§_= l/a — s- 
0), the nearest-neighbour operator Hi in (1.1) acts as a 
perturbation operator and we therefore have to compare 
in first order t he e xpect atio n values of Hi between the 
ground states (A2) and (A7), respectively: 
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{^{s,ps = (0,0),iV)|ifi|*(s,p, = (0,0), iV)) 

= J2 (*<+H<T)|53(x)|vI'(+)(a)) 

xe(+) 

X ^(vI/(-)(a)|T,53(x)r-i|vI/(-)(a)) . (AlO) 

Here we have used the fact that the nearest-neighbour 
couphngs in Ti\ connect the spin operator S(x) on the 
even sub-lattice with the spin operators S(y = x -I- z/) = 
Tjj^{x)T~^ on the odd sub-lattice. 

According to (A3) and ( [A5| ) the application of the 
translation operator yields: 



^(*(-)(a)|r.53(x)T-i|vI/(+)> 
4- (*W(a)|53(x)|vl/W(r7)) 



(All) 



7V/2 



■ 



G = s/2 , 



= 4M. 



We end up with the following expression for the ground 
state energies in the antiferromagnetic collinear phase 
a > ao{M) in the first order approximation: 

E{s, p, = (0, 0), a, N) = a[2E{s/2, p,/2 = (0, 0), 7V/2) 



1 



2Af2 • N 



(A12) 



We compute in the same manner the exp ectation val- 
ues of Hi between the ground state ( |A7| ) for s odd, 
a = s/2 + 1 even (p^ = (0, tt)): 



E{s, Ps = (0, tt), a, N) = a[^E{a, p, = (0, 0), 7V/2) 

+E{a -l,p, = (7r,^),iV/2) 



+ -2M^N 
a 



(A13) 



APPENDIX B: SECOND ORDER 
PERTURBATION THEORY IN THE STRONG 
FRUSTRATION LIMIT 

The evolution of the energy eigenvalues and transition 
matrix elements of a Hamiltonian 



i?2 + f3Hi 



(Bl) 



depending on some parameter j^is described by a closed 
system of differential equation^^. In particular one finds 
for the second derivative of the ground state energy: 



fEo 
df3^ 



En — Eq 



(B2) 



T„o = (n|iJi|0) 



(B3) 



|0) and In) denote the ground state and the exited 
states of (Bl) with energy eigenvalues Eq and _E„, re- 
spectively. 



We are now going to evaluate (B2) in the strong frus- 
tration limit /3 = 1/a ^ oo. The ground state of H2 in 
the sector with total spin s (s/2 even) is given in (|A2|). 
In the following we will consider those contributions to 
the transition amplitudes, which arise from the excited 
states: 

**(s,p, = (0,0),A) = 

(a + 1, p, N/2)¥-\a - 1, -p, A^/2) 

+ *(+)(ct - 1, -p, A/2)1'(")(cr + 1, p, A^/2)) (B4) 

with energy i;(+)(cr + l,p, A^/2) + £;(~)(cr - l,p, A^/2) 
and a ^ s/2. E'^^^a' ,p, N/2) are the lowest energies 
on the subsystems for a given spin a' and momentum p. 
The momenta p and — p on the subsystems have to add 
up to zero in order to produce a translationally invariant 
(momentum Ps = 0) state on the whole lattice. Similarly 
the sublattice spins a + I and a — 1 add up to the total 
spin s of the whole system. Any other distribution of 
the total spin s on the sublattices will increase the total 
energy in comparison with i?^+^((T+l, p, N/2)+E^~'>{a— 
l,p,N/2). 

The computation of the transition matrix elements 
yields: 



= E (<\.pl^+(^)l*il=(o,o)) 
xe(+) 



--{0.0)' 



o-,p„ = (0,0)/ 



N 
'2 



x(<\,p|5+(x = 0)|vI/W^(„„)) 
<(M^W^_p|^-(x = 0)|vI/W^(„_„^) 



(B5) 



Here we have used the transformation properties ( A5) 
of the sublattice eigenstates 4''^+)((t ± l,p) under the 
translation operators Ty. 

A similar expression for the transition matrix el em ents 
to the second component in the excited state (B4), is 
obtained by the substit utio n p —* p in the prefactor on 
the right-hand side of (B5). 

Combining both results, we obt ain for the transition 
probabilities to the excited states (B4): 



where 
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i(*:,pj^imF 



(0,0) 



f N\~ 2 2 

2(yj -(l + COSpi) (I + COSP2) 

x|(vI/W,p|5+(0)|*(^^^(„^„))P 



For the evaluation of the second order term (B2) in the 
perturbation expansion of (Bl) around (3 — l/a = we 
also need the energy differences 

En-Eo^ En, (P, (7 + 1, iV/2) + En, (-p, (T - 1, iV/2) 

-2Eo{po^O,a,N/2). (B7) 

We have computed the energy differences 

En, (p, fT ± 1, iV/2) - Eo{po = 0, a, N/2) , 

a = s/2, 1^1,2 

on the zb-subsystems as well as the transition proba- 
bilities 

l(*itl±i,±pl^±(0)l*-,p.=o)P (B8) 

by means of the recursion mcthotS. The resu ltin g 
coefficient (5_i(M), which appears in the expansion (2.4): 



is shown in Fig. pXl. 



(B9) 
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